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ABSTRACT 


Three double-gimbnled control moment gyros constitute the actuators 
for one of the altitude control system modes of the Apollo Telescope Mount. 

A distribution law and a rotation law for the individual angular momentum 
vectors of the control moment gyros arc developed because the attitude control 
law, in controlling the magnitude and the direction of the total angular momentum 
vector, uses only three of the available six degrees of freedom. Without a 
distribution law, a highly undesirable distribution can develop where two of 
the three individual angular momentum vectors are parallel and the third is 
untiparailel, reducing the total angular momentum available to one third of 
that which could be available with a desirable distribution. The desirable 
distribution is such that, after equilibrium is achieved, the individual angular 
momentum vectors have equal components along the total. Other results are 
equal angles between the individual vectors and equal angles between the 
individual vectors and the total. Because of the equal angles, the desired 
distribution is called the "isogonal distribution. " Isogonal distribution has 
the following advantages: (1) antiparallel distribution is avoided, (2) the gain 
available for the attitude control law is maximized, and (3) the cross coupling 
is minimized. The distribution law uses two degrees of freedom; rotation 
about the total angular momentum vector still remains free. Therefore, a 
rotation law is added to minimize ihe inner gimbal angles with the result that 
hitting of the inner gimbal stops is avoided ns much as possible. Operation 
of the distribution and rotation laws in case one of t he control moment gyros 
fails is also discussed. 
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DEFINITION OF SYMBOLS 


i, J, k 


m, n 


B. 


cyclic permuted subscripts (1,2,3, or 2,3, 1 or 
3,1,2) 

general subscripts (range 1,2,3; all possible 
combinations) 

a bar below a letter indicates a vector quantity 
a dot above a letter indicates a time derivative 
a 2x3 matrix (eq. 9) 
a 3x2 matrix (eq. 11 and 11a) 


- cos 


% 


a unit yector in direction of the ith CMG angular 
momentum 


mn 


direction cosines of e ^ (i = m) with respect to the 
vehicle axis n 


'ik 


-T 


= n / 1/(1 - e . k 2 ) 

“ e t + e^ + e 3 . Normalized total angular momentum 
of the CMG cluster ( maximum value: 3) 


= I e T | 


Ti 


components of e T in vehicle space 


intermediate quantity used in the digital program 


iv 



DEFINITION OF SYMBOLS (Continued) 




mn 


K. 


k d' 


K 


R 


mn 


0 -e 3 +©j 

+e 3 0 - ei 

-e 3 +Cj 0 


see eq. 5 

components of (i = m) in vehicle space 
gain of the isogonal distribution, [ 1/s] 
maximum of IK^I , ( l/s] 
rotation gain (K^ > 0) , [ 1/s] 
determinant (see eq. 13) 

intermediate quantities used in the digital program 


= sin 


t 


*i 



a. 


I 

\ 


- tan 

unit vector along the ith vehicle axis 
vehicle axes 

angles used in conjunction with the rotation law, 
defined in Figure 4 
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DEFINITION OF SYMBOLS (Concluded) 


inner gimbal angle 


outer gimbal angle 
6 l(i) 

L a 3(i>_ 


of ith CMG 
of ith CMG 


error of the isogonal distribution, [ 1/s] 
rotation error, [ 1/s] 

part of e D stemming from the ith CMG, [ 1/s] 

ri 

= e T 2 

gain modifier for isogonal distribution 
actual angular velocity of e ,, , ( 1/s] 
angular velocity command for rotation, ( 1/s] 
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INDIVIDUAL ANGULAR MOMENTUM VECTOR DISTRIBUTION 
AND ROTATION LAWS FOR THREE DOUBLE-G IMBALED 
CONTROL MOMENT GYROS 


SUMMARY 


A distribution law and a rotation law for the individual angular momentum 
vectors of the three control moment gyros of the Apollo Telescope Mount are 
developed. The attitude control law controls the magnitude and the direction of 
the total angular momentum vector, using only three of the available six degrees 
of freedom. Without a distribution law, a highly undesirable distribution can 
develop where two of the three individual angular momentum vectors are parallel 
and the third is antiparallel, reducing the total angular momentum available 
by two thirds. The total is always fully available with the desirable distribution 
as provided by the distribution law, unless gimbal stops interfere. ihe desirable 
distribution has the following advantages: ( 1) antiparallel distribution is 
avoided, (2) the gain available for the attitude control law is maximized, and 
(3) the cross coupling is minimized. The rotation law minimizes the inner 
gimbal angles with the result that hitting of the inner gimbal stops is avoided 
as much as possible. Operation of the distribution and rotation laws in case 
one of the CMG's fails is also discussed. 


INTRODUCTION 


The Apollo Telescope Mount (ATM) carries three double -gimbaled con- 
trol moment gyros (CMG) [ 1,2] . The magnitude of the angular momentum of 
each CMG is fixed, and the arrangement has six degrees of freedom. The 
attitude control law will command the magnitude and the direction of the total 
angular momentum; therefore, it utilizes only three degrees of freedom. The 
distribution of the individual angular momenta is not controlled. This freedom 
can lead to highly undesirable distributions, such as the antiparallel situation 
in which two vectors are parallel and the third is antiparallel. To avoid un- 
desirable distributions, a distribution law is developed. The desirable 



distribution is such that, after equilibrium is achieved, the individual angular 
momentum vectors have equal components along the total. Other results are 
equal angles between the individual vectors and equal angles between the 
individual vectors and the total. Because of these equal angles, the desired 
distribution is called the "isogonal distribution" (Fig. 1). The isogonal dis- 
tribution has the following advantages: ( 1) antiparallel distribution is avoided, 
(2) the gain available for the attitude control law is maximized, and (3) the 
cross coupling is minimized. The distribution law uses two degrees of freedom; 
rotation about the total angular momentum vector still remains free. There- 
fore, a rotation law is added to minimize the inner gimbal angles with the 
result that hitting of the inner gimbal stops is avoided as much as possible. 

The existing gimbal angles are the inputs to the distribution and rotation 
(D&R) laws and the gimbal angle rate commands are the outputs. These rate 
commands are in addition to those from the attitude control law. To develop 



FIGURE 1. ISOGONAL DISTRIBUTION 
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the D&R laws it is assumed that the angular velocity of the vehicle is negligibly 
small and that the contribution of the gimbal rates caused by the D&R laws 
to the total angular momentum is negligible. For simplicity, unit vectors 
along the individual angular momentum vectors are used and are called 
e. (i = 1, 2, 3); their sum is e T (0 t |e^,| ^ 3). To avoid interference of 

the D&R laws with the attitude control law, it is necessary that these laws do 
not result in any change of the total angular momentum vector of the CMG 
cluster (i.e. , no acceleration on the vehicle). 

The ratio of the angular velocity command to an angular error function 
may be called the gain for the D&R laws. In addition to a selectable constant, 
there is an inherent gain which varies with the gimbal angles. This is accept- 
able if the loop is stable with the highest inherent gain. 


DISTRIBUTION LAW 


For the sake of argument we can assume that two of the e ^-vectors are 
summed first and then the third is added to form e If we desire to redis- 
tribute the first two without disturbing the total (e = 0) , we can only rotate 
the first two about their sum. The angular rate e about this sum is made 
proportional to the scalar function e. which, for e t and e 2 , is given by 
Kple j . e T - e 2 . SLj)* This scalar function provides a signal proportional 

to the deviation from the desired orientation. To assure rotation about tne sum 
of the first two vectors, e is also made proportional to their cross product. 

Pairing ej and e 2 , for example, we have (Fig. 2) 

£i = £3 (e 2 v e,) ( 1 ) 

e 2 = e 3 (e j x e 2 ) (2) 

with e 3 = K d (e t • e T - e 2 * £ T ) = K^e, • e 3 - e 2 • e 3 ) (3) 

By cyclic permutation the equivalent equations for the other two pairs 
can be developed. The constant can be chosen from other considerations 

(see discussion at the end of this section) . The cyclic permutation allowr, us 
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FIGURE 2. DEVELOPMENT OF DISTRIBUTION LAW 
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to use index notation with the understanding that we have the following three 
sets of values for i, j, and k: 


i 

i 

k 

1 

2 

J 

3 

l 

2 

3 

1 

3 

1 

2 


The vector pairing can be assumed tc be simultaneous and the resulting 
velocities vectorially added to form: 


e . = h . x e . 

-i -i -i 

(4) 

h . = e, e. + e.e, 
-i k-j j-k 

(5) 


It can be seen ( after expansion of equations 4 and 5) that no matter 
what value is chosen for K^, the sum of the e^'s (ore^,) is zero (Appendix A). 

To develop the equations for the gimbal velocity commands, we must 
consider the following built-in relationship (Fig. 3): 


e.. 

li 


e . = 

— l 



cd,. .c<5„..,u. - c5.. ..s6„. .. u . - s<5.,.. u 
l(i) 3(i)~i l(i) 3(i)~j l(i) — k 


( 6 ) 


where the u ,'s are unit vectors along the vehicle axes with the indices permuted 

as previously shown, and the 6's are gimbal angles defined in Figure 3. The 
second index on the components of e . refers to the vehicle axes; for example, 

e. is the direction cosine of e . with respect to vehicle axis j. Therefore we 
i] -i 

can also write: 
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e ll e l2 e 13 M.1 

e 21 e 22 e 23 H.2 

e si e 32 U 3 


The generally nonorthogonal direction cosine matrix in equation 7 will be 
needed iater. The e„'s are gi\en in terms of gimbal angles on page 20. 

For an existing gimbal velocity we have 


e. = A .6. 

-i i— (i) 


e,. and 6... = 

ij “(l) 


Ad 


1 ( i) C °3(i) 

-c6 w .. s<5„, .. 
l(i) 3(i) 

'l(i) S<5 3(i> 

-c6 , c6 . . 
1 ( i) 3(i) 

l(i) 

0 


As 6 .-command law we need (see Appendix B and proof below) : 


6,., = B. h. 

^i) i-i 


+S6 3U) 

+c<5 

■ t6 l(i) Cd 3(i) 

+ “i 
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or in terms of direction coeines (with e^' 


= n / 1/(1 


ik 


) ): 



" e ik' e ij 
e ik (e ik )2e ii 


e 'e 
ik ii 


- e (e *) 2 e 
ik ik ij 



( ila) 


The existing hardware implementation through resolvers on the gimbals provides 
the direction cosines of equation 7 whereas the gimbal angles are not available 
explicitly. Thus, IT is expressed in terms of the direction cosines (see 
Appendix C). 

By multiplying, it can be shown that 


e = A. B. h . = -e. xh=h,xe, 

-i ii-i — i — i — i -i 

if it is assumed that the actual and the commanded gimbal rates are identical. 
Equation 10 is, therefore, the proper solution for the gimbal velocity commands 
to solve equation 4. 

A discussion on the freely selectable isogonal distribution gain (equa- 
tion 3) is in order. When this gain is positive, we get a right-handed con- 
figuration ( right-handed means that when we look down e ^ and identify the 

individual vectors counterclockwise, we find the sequence 1-2-3); for a 
negative K^, we get a left-handed configuration. Both configurations are stable. 

Since neither sign has a basic preference, we can dynamically (during operation 
of the CMG's) select the sign of the gain so that the vectors are driven 

toward the closer solution. This action implies a knowledge of whether the 
instantaneous distribution is right- or left-handed. The following determinant 
affords this knowledge. 


e u 

e 12 

e l3 

e 21 

e 22 

e 23 

e 31 

e 32 

e 33 
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It can be recognized that this is the determinant of the direction cosines for 
the angular momentum vectors. If these vectors are mutually perpendicular, 
we nave an orthogonal coordinate system, and it is well known that the value 
of the determinant q is +1 for a right-handed system and -1 for a left-handed 
one. For a continuously deformable system as the one at hand, the detr 
minant can be anywhere between these limits and the sign of the determinant 
can be used as the sign for the distribution gain K^. A value of zero for q 

indicates that the three vectors are in a plane. A change in sign for q occurs 
under certain conditions. Gimbal stops may disallow the isogonal distribution 
and deform the distribution so much that q changes sign. The distribution gain 
is much lower than the attitude control gain (as discussed later), and maneuver 
commands can cause a relatively fast change in the distribution which could 
resuit in a sign change of q. Finally, whenever e ^ passes through zero, q 

changes sign. Regardless of the reason for the sign change, it is always bene- 
ficial for the system behavior. A further improvement in system behavior 
can be made through the introduction of a multiplier A^ so that the isogonal 
distribution gain becomes 


k d = *d k d' sgn 1 

where 

K_' = | K | 

l> D max 

and 

0 S } < +1 

The multiplier A^ is a function of e^ = |e ^1. For small e^ we desire A^ = 0, 

thus avoiding the possibility for large angular redistributions when e ^ passes 

close to zero. When e^, is in the vicinity of unity (i.e. , without a distribution 

law, antiparallel is a problem), the multiplier A^ should be at its maximum. 

For large e T , a small A^ is sufficient. The exact function of ^ D ^ e T ) is not 

critical, and straight line segments may be used (Appendix C). The maximum 
distribution gain K^' should be large enough to avoid antiparallel, but not 

larger than necessary because gimbal stops sometimes disallow the isogonal 
distribution and the attitude control gain must be higher than the distribution 


9 



gain so that the attitude control takes over in case of conflict. During normal 
operation there is no conflict because the distribution and rotation laws do not 
result in an acceleration of the vehicle. 

Temporary deviations from the isogonal distribution can be tolerated to 
a large extent. For example, for the antiparallel case, e^, = 1. This case is 

depicted in Figure 1. The dot product between any of the isogonally distributed 
vectors and the antiparallel direction is - 1/3, or there is an angle of almost 
110 degrees between the desired individual vector direction as demanded by 
the isogonal distribution law and the antiparaliel direction. This indicates 
that even a deviation of 40 or 50 degrees from the desired direction can be 
temporarily accepted. 


ROTATION LAW 


The isogonal distribution law concerns itself only with the relative 
distribution between the individual and the total angular momentum vectors. 
A rotation about the total leaves the distribution unaffected. The rotational 
freedom can therefore be used for some benefit. To reach some desired 
orientation about e^, an angular velocity of the following form must be 
commanded. 

0 = e R e T (14) 


where e will be developed later. For the individual angular momentum 
K 

vectors, the relation exists 

e. = !]xe. (15) 

Equation 15 is identical in form to equation 4; only Q is substituted for h.. 
Therefore we get analogous to equation 10 



(16) 


No unique law exists for the rotation contrary to the unique isogonal 
distribution law. In the application of the three CMG system to ATM, the 
most benefit can be derived from a rotation law if the inner gimbal angles 
are minimized. Even under this assumption, minimizing must be first defined, 
since no unique definition exists. Avoidance of hitting the inner gimbal stops, 
which are at ±80 degrees, was given priority. It was also desirable to have a 
continuous function rather than a switching function for the rotation law. Be- 
cause the three vectors rotate as a unit about e^, none of the inner gimbal angles 
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will generally be at an absolute minimum, but a compromise will be reached 
for all three. The following form was selected for : 


f R = 


3 

Z e 


i=l 


Hi 


with 


k r 

e Ri = 7 ^ 


(17) 


( 18) 


Because will become zero at equilibrium and K is a freely selectable 
positive constant, no additional constants a.re necessary in equation 17. 


The reasons for the various quantities in the error function are 

given. To make ft independent of the magnitude of e ^ we have lo divide 

by e^. The tangent of the inner gimbal angle ( rather than the angle itself) is 

chosen for two reasons: (1) to give the larger angles more weight and (2) to 
make use of trigonometric functions which are readily available in the form of the 
direction cosines (equation 7). Figure 4 defines the two other angles used for 
i = 1. Consider, for the sake of the explanation, that the plane for which 
<5 1 ( i) = 0 is the equatorial plane. Then the angle a. indicates the separation 

of the meridian plane containing e . and the one containing e^. The combination 
of tfi^^sa. will select the proper sign for the angular velocity and go to zero 
when the desired conditions are reached. The multiplier c/3, takes the effective- 
ness of a rotation in reducing 6 . , .. into account, eliminating the e -contribution 

l(i) Ki 

to completely when 0. = where a rotation does not change . 

We can identify several cases; Figure 5 shows an opaque unit sphere, the 
equator of which indicates 6^ = 0. The vertical line is a longitude onto which, 

for convenience, all intersects for the direction of e T have been placed. Usually 

e^, will not terminate on the unit sphere. 


CASE I . - 
rotation about e rf ) 
<5^^ is positive, 
unstable. 


The circle of possible tip locations ler e . (assuming a 
does not intersect the equator, and the inner gimbal angle 
We find two equilibrium conditions: one stable and one 
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S|(i)-+90* 



TL» CIRCLE OF U = UNSTABLE 

TIP LOCATIONS OF ej 8, (|) --90° EQUILIBRIUM 

FOR ROTATION S = STABLE 


ABOUT e T EQUIUBRIUM 
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Determination of equation 18 in terms of the readily available di. ection 
cosines and the components of the total angular momentum vector leads to 
(i = 1). 



e n e T2 " e i2 e Tl 
n /(1 - e" l3 2 ) (e T i 2 + e T 3 


e Rl = "c” 2 j" e t 3(ei3 ' )2 (eil6 T2 ” ei2G Tl ) [ 

Cyclic permutation of the indices leads to e R2 and e R3 : 

e R2 = “ e 21 ^ e 2l') 2 ^ e 22 e x3 ~ e 23 e T2*| 

€ R3 = e^ 2 ) “ e32 ^ e3: ^ 2 (e 33 e Tl ” e3lG T3\' 


( 19 ) 


( 20 ) 


( 21 ) 
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DISTRIBUTION AND ROTATION LAW COMBINATION 


Since the desired gimbal velocity commands from the distribution law 
and from the rotation law are independent of each other, they can be combined 
into 

4(i) = W VLt' (22 ' 

These gimbal velocity commands should be added to the commands from 
the attitude control law. 

The distribution is unaffected by the rotation as already indicated during 
the development of the rotation law; i.e. , the e.'s of the distribution law are 

not affected. On the other hand, e is affected by a change in tl.s distribution, 
and it is advisable that (contained in e^) is made smaller than K^’ {con- 
tained in h. through e.). A ratio of K /K ' = 0.1 was found to be acceptable, 
l i K D 

but it is not critical either way. 


OPERATION WITH TWO CMG'S 


A system with three CMG's has the capability to lose one and still be 
able to control the vehicle, even if the performance is degraded. Two CMG’s 
have only four degrees of freedom, and three are being used for attitude con- 
trol. There is no need for a distribution law since the two remaining angular 
momentum vectors have inherently the proper distribution, but the rotation 
law still can minimize the inner gimbal angles of the operative CMG's. To 
convert from the three to the two CMG operation, it is only necessary to set 
all the direction cosines of the inoperative CMG to zero. (Gimbal angle velocity 
commands issued by the computer to the inoperative CMG should be dis- 
regarded. ) No distribution command will appear for the operative CMG's. 

But, as desired, the rotation law is working as usual, having now only two 
error sources, rather than three (equation 17) . 
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CONCLUSIONS 


The isogonal distribution lav/ and the rotation law were implemented in 
a hybrid simulation, and the behavior was studied with the aid of a three- 
dimensional display. ♦ Both laws performed as expected. It was also noted 
that, if gimbal stops disallowed the isogonal distribution, the distribution 
approximated the desired distribution as close as possible under the imposed 
restrictions. Analysis cf the dynamic performance of the attitude control law 
showed a marked improvement in response of the desired channel and a strong 
reduction of the cross coupling into the other channels. 


♦The three-dimensional display shews the body axes, the three individual 
angular momentum vectors, their total, and the total commanded angular 
momentum vector simultaneously, identified by appropriate coding. This 
proved an invaluable tool for the investigation of the behavior of the system f 3] . 
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APPENDIX A 


Sum of the e's 

Expansion and combination of equations 4 and 5 yield 
e,= (eje, + e 2 e 3 ) x e t 

ki = (€1^3 + € 3 e,) x e 2 (A-l) 

e 3 = (€ 2 ej + ) x e 3 

For the sum of the e/s, we get 

e T = + (§2 x e 3 + e 3 x e^ ) 6! 

+ (e 3 x ej + e 4 x e 3 ) e 2 (A- 2 ) 

+ (ij x §2 + e, x e t ) €3 

Because the vectors in parenthesis are all zero, e^ is zero and it is immaterial 
what is chosen for the s.'s. The fo^m chosen was (equation 3 ): 

€1 = K D (e 2 - e 3 } • e* 

€ 2 = K D ^3 ' > ' % 

€3 = K D (ei - §2 ) • e 3 

Because the s.'s are indenendent of each other, the K 's could be different for 
1 D 

each t. without affecting the condition that e_ = 0 , This is not sensible though, 

1 0 -T 

because it would no Teat the vector pairs alike, and there is no reason for that. 
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APPENDIX B 


Development of Matrix B 


We desire to get the gimbai velocities to solve equation 4, which is 
(subscripts dropped) 


e = h x e 


(B-l) 


The actual angular velocity of e can be defined as uj and must be per 
pendicular to e because of the gimbai constraints. Then, 


e = x e 

This leads to (equations B-l and B-2) 


(B-2) 


(h - a;) x e =0 

and this equation can only be solved because cu and e are known in terms of the 
gimbai angles and their velocities. With 



" e i“ 


+c5 1 c<5 3 


+s6 3 5, 

e = 

e 2 

= 

-c6jS6 3 

and UJ = 

+ 063 ^ 


_ e 3_ 


-s 6 j 


_ -«l - 


we obtain for equation B-3 



0 -(h 3 + <5 3 ) + (h 2 - c6 3 5,) 


+c6 s c6 3 

0 = 

+ (h 3 + 6 3 ) 0 -(hj - s6 z 6i) 


-C 6 JS 65 


-(h 2 - c6 3 6 2 ) +(hj - so 3 6j) 0 


- ~s6t _ 
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The third equation in 3-4 yields 


6j - +s63hj + c<5 3 h 2 ( B-5) 

With equation B-5 the first or second equation of B-4 can be solved to get 

63 = +t<5 1 (s6 3 h 2 - c6 3 h,) - h 3 (B-6) 


or 



It might be interesting to note that the cross-product law is 

(XPR) = BAB 
and that 

e = - A B 

where e is defined as 
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APPENDIX C 


Suggested Digital Computer Implementation 
(Sequence of Operations) 


Inputs. 


e ll " 

+Ca i(l) Ca 3U) 

k d' 

e l2 = 

■ C6 1(1) S6 3(1) 

Kr 

e 13 = 

- S4 1(1) 


e 21 = 

- S6 1(2) 


e 22 = 

+Ca K2)° 6 3(2) 


e 23 = 

' C6 1(2) S6 3(2) 


e 31 = 

_c6 1(3) s6 3(3) 


e 32 = 

- S6 1<3> 


e 33 = 

+c6 1(3)° 6 3(3) 


The equivalences of the e 1 

mn 

s in terms of the gimbal 


given for completeness. The e 's are provided by resolver chains in the 
existing hardware and constitute the actual inputs to the digital computer. 
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Sequence of Computations . 


(e i3 ') 2 = (l-eu 2 )* 1 1 

both quantities are 

e„' - l<e, 3 ') 2 ] +i \ 

needed in the following 

<e 31 ') s = U - e,, 2 )-' 1 


e n ' = [ (e 31 ') 2 J + 2 \ 


( e 32*) 2 = ( f ~ e 32 2 ) 1 I 


\ 

> 

e s! ' = [<e 3! ’) 2 i + 2 \ 



e Tl = 6n + 021 + 631 

e T2 = e 12 + e 22 + e 32 

0 T3 = 013 + 623 + 633 


p 2 = « 2 0 2 + g 2 

T T1 T2 T3 


e T ' / *7 


q = e 11 (e 2 2e n - 632623) + e2i(e 3 2e is - e 12 e33) + e 31 (e 1 2e23 - 6326!,) 
sgn q = + i for q = 0 

sgn q = -1 for q < 0 
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= 0 for 5 0. 25 

A d = 2e T - 0.5 for 0.25 < = °- 7 5 

A d = +1 for 0.75 < e T S 1.25 

A d = +3.5 - 2e T for 1.25 < = 1-65 

= +0.2 for 1.65 < 

k d = k d' *d sgn « 

E 1 = e 21 e 31 + e 22 e 32 e 23 e 33 

E 2 = e u e 31 + e 12 e 32 + e 13 e 33 

E 3 = e ll e 21 + e 12 e 22 + e 13 e 23 

Cl = (E 3 - E 2 )K d 
e 2 = (El - E 3 )K d 

e 3 = ( E 2 “ e i)Kjj 

h ll = e 3 e 21 + e 2 e 31 
^12 = e 3 e 22 + e 2 e 32 
hj 3 = e 3 e 23 + e 2 e 33 

h 21 = e l e 31 + e 3 e ll 
^22 = e l e 32 + e 3 e 12 
^23 = e l e 33 + e 3 e 13 


22 



h 3i = € 2 e li + € i e 21 

^32 = ^2 e 12 + e l e 22 

^33 = e 2 e 13 + € l e 23 


€ R1 I = e 13 (e 13 f ) 2 (e la e T1 - e^e^) 

= e 2i^ e 2i r ) e 23 e x 2 * e 22 e T3 ) 


^g' “ e 32^ e 32') 2 ^ e 3i e T 3 “ e 33 e T1 ) 


» 8 

T 


e R ~ e J * € R1* + W + W* 


I'll = h ll + 

€ R e Tl 

r 12 = 1 »,2 + 

€ R 6 T2 

*13 = h 13 + 

C R e T3 

1 

to 

►* 

II 

jsr 

+ 

e R e Tl 

r 22 = h 22 + 

€ e 
R T2 

r 23 = h 23 + 

e R®T3 

*31 = h 31 + 

e R®Ti 

r 32 “ h 32 + 

e R®T2 

r 33 = h 33 + 

e R®T3 
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Outputs. 


<5^^ = - e 13^ e 12 r ll ■ e ll r 12^ 

= + ( e 13') 2e 13^ e ll r ll ‘ e 12 r 12) ~ r 13 

< 5 ^ 2 j = -e 2 i'(e 2 5 r 2 2 - <2221*23) 

^3(2) = + ^ e 2i') ?e 21^ e 22 r 22 + e 23 r 23 •) ” r 21 

= ”632' ( e 3l r 33 " e 33 r 3l) 

63(3) = + (es 2 ') 2 e 32 (e33r33 + e 31 r 31 ) - r 32 
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